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Abstract 

Period doubling bifurcation leading to subharmonic oscillations are undesired phenomena in switching 
converters. In past studies, their prediction has been mainly tackled by explicitly deriving a discrete time 
model and then Unearizing it in the vicinity of the operating point. However, the results obtained from 
such an approach cannot be applied for design purpose. Alternatively, in this paper, the subharmonic 
oscillations in voltage mode controlled DC-DC buck converters are predicted by using a formal sym- 
bolic approach. This approach is based on expressing the subharmonic oscillation conditions in the 
frequency domain and then converting the results to generalized hypergeometric functions. The obtained 
expressions depend explicitly on the system parameters and the operating duty cycle making the results 
directly applicable for design purpose. Under certain practical conditions concerning these parameters, 
the hypergeometric functions can be approximated by polylogarithm and standard functions. The new 
approach is demonstrated using an example of voltage-mode-controlled buck converters. It is found that 
the stability of the converter is strongly dependent upon a polynomial function of the duty cycle. 
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I. Nomenclature 

Let us define the following parameters related to the buck switching regulator studied in this work and 
that are listed in analphabetic order. 

^ _ '^^Jz2 Ratio between ll)z2 and ojq 

C Output Capacitance 

D Steady state duty cycle 

fs Switching frequency 

Uzi = --^ Zero due to the ESR rc 

kv Proportional gain 

L Inductance 

= C{Rr. + Rr, + r.n) + L ^^'"P^"^ ^^^^^^^^ '''''''' 
R Load resistance 

Tc ESR of the capacitance C 

ESR of the inductance L 

T Switching period 

Vg Input source voltage 

Wo = \ / 7TT— — Angular resonance frequency 

V {R + rc)LC 

LOip Pole of the type II controller 

LOs = 27r/s Angular switching frequency 

ujzi Zero due to the ESR of output the capacitor C 

L0z2 Zero of the PI controller 

II. Introduction 

Switching converters are integral elements to modern power electronics. Despite their widespread 
use, they can pose serious challenges to power-supply designers because almost all of the rules of thumb 
governing their design are only applicable to the linearized averaged system even though the system works 
in switched mode |[T|, Q. Switch mode operation is carried out by means of Pulse width Modulation 
(PWM) action on the switches. In the traditional PWM control, the duty cycle of the pulse driving signal 
5{t) is varied according to the error Ve{t) between the output voltage and its desired reference. This 
error is processed through a compensator to provide the control voltage Vc{t). The simplest analog form 
of generating a fixed frequency PWM is by comparing the control voltage with a ramp periodic signal 
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(a) 



(b) 



Fig. 1. (a) Block circuit diagram a DC-DC buck converter under voltage mode control with a controller with transfer function 
Gc(s) and a PWM modulator, (b) Equivalent block diagram. 



vtri{t) in such a way that the pulse signal goes high/low when the control signal is higher/lower than 
the triangular signal. A rising ramp carrier will generate a PWM signal with trailing edge modulation 
and a falling ramp carrier will generate a PWM signal with leading edge modulation |3|. In the trailing 
edge modulation strategy, when the control voltage is greater (resp. smaller) than the sawtooth ramp 
voltage, the switch is ON (resp. OFF). Whereas in the leading edge modulation strategy, the opposite 
happens. Figure [T] shows a schematic circuit diagram of a buck converter under voltage mode control 
with a compensator Gc{s) and a fixed frequency PWM and its equivalent block diagram. 

Probably, the buck converter is the most common voltage regulator in use. Its structure is not com- 
plicated (See Fig. [T]-(b)). The power stage together with the modulator consists of a pulse generator and 
a passive low pass filter. For this reason, it can be used to convert an input source voltage Vg into a 
lower output voltage v. In order to regulate the output voltage, a voltage reference, an error detector, a 
compensator and a modulator are added to the circuit. These basic elements form a complete switching 
buck DC-DC regulator. 

In practice, it is desirable that the system operates periodically with a constant switching frequency 
fs = 1/T equal to that of the external sawtooth ramp modulating signal which is the same frequency of 
the external clock signal. However, under parameter changes, the stability of this operating mode may 
be lost resulting in different kinds of instabilities and dynamical behaviors Q- |12|. There have been 
hitherto many research efforts devoted to predict the border of occurrence of such instabilities. However, 
in the case of subharmonic oscillations in voltage mode control, most of the reported analysis are based 
on abstract mathematical models, which are demonstrated to be very powerful to accurately predict this 
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phenomenon, but not so suited to obtain clear design-oriented criteria in the parametric design space. In 
the past studies, the prediction of subharmonics is mainly based on deriving an accurate discrete time 
model and then linearizing it in the vicinity of the operating point. However, the results obtained from 
such an approach cannot be applied for design purpose. Recently, Filippov's method and the monodromy 
matrix were used to predict these instabilities in DC-DC converters and similar results to those obtained 



from the discrete time approach were derived |13|, |14|. 

While there are powerful design-oriented techniques to characterize low frequency instabilities related 
to the averaged model ||T|, there is still a lack of design oriented tools for predicting fast scale instability 
in the form of period doubling and subharmonic oscillations inherent to the switched nature of the system. 
Note that this phenomenon can be perfectly predicted in the peak and the valley current-mode-controlled 
system examining mainly the slopes of the inductor current ( yj) or by using modified averaged models 
taking into account delay terms due to the sample-and-hold effect ( p3) ) but not in voltage mode control 



or in average current mode control |16|, |17|. For example, it is known that to guarantee the stability 
of a peak current-mode-controlled buck converter, the following inequality, in terms of the duty cycle D 
and the system parameters, must hold 

«-i<^ (!) 

2 Vg 

On the other hand, voltage mode controlled power supply has become very popular especially in low 
noise next generation communication systems. Such systems are supplied by means of different Point of 
Load (PoL) voltage regulators in the form of single phase voltage-mode-controlled buck converters |18|- 



1 19 1. It is well known that these converters, under voltage mode control, are prone to exhibit undesired 
subharmonic oscillations for certain parameter values Q, |[6|, ||8|. Figure |2] shows typical waveforms and 
state plane trajectories for the desired periodic behavior and for subharmonic oscillations in a voltage- 
mode-controlled buck converter. 

An equivalent expression to ([T]), in the case of voltage mode control, will be extremely helpful in 
designing switching converters free from subhamonic oscillations. One may, therefore, ask: There exists 
a similar expression as ([T]) for the case of the voltage-mode-controlled system? This paper tries to develop 
such an expression for different voltage mode control schemes. 

The rest of this paper is organized as follows: a short tutorial on hypergeometric series and their 
special case, polylogarithm functions, is given in Sections III. Section IV will revisit the dynamic model 
of the buck converter in the frequency domain and conditions for periodic behavior and subharmonic 
oscillations are presented in this domain by using Fourier series. These series are calculated exactly by 
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(e) (f) 

Fig. 2. Different dynamical behavior of a DC-DC buck converter under voltage mode control. Desired behavior (left) and 
subharmonic oscillations (right), (a)-(b) Time domain waveforms of the control voltage Vc(t) and the ramp voltage vtri{t), 
(c)-(d) corresponding state plane trajectories, (e)-(f) corresponding experimental data. 
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using hypergeometric functions and approximately by using polylogarithmic functions. The approach is 
then appUed, in Sections V, to a voltage-mode-controlled DC-DC buck converter with different controller 
transfer functions. Some circuit-based simulations are presented to validate the results obtained from the 
derived theoretical results. In Section VI, the results are formulated in terms of Figures of Merit, widely 
used by the power electronic community, such as the crossover frequency and the phase margin. Finally, 
in the last section, some concluding remarks of this work are summarized. 

III. A Short Tutorial on the hypergeometric and polylogarithm functions 
A. Gamma function 

It is somewhat problematic that a large number of definitions have been given for the gamma (F) 



function 1 20 1 . Although they describe the same function, it is not entirely straightforward to prove their 



equivalence. One of the definitions of the gamma function due to Euler is |20| 



F (z) = lim 



1 -i-r 1 + n 



'-U'-^ (2) 



n-5>oo z(z + 1) . . . (z + n) z z + n 

n=l 

Using integration by parts, one can show the recursive relation T {z + 1) = zF (z) in such a way that for 
integer numbers n, F (n + 1) = nl. The Euler gamma function is one of the most important functions 
in mathematical physics pO} . It may be regarded as a generalization of the factorial which spreads its 
values over the whole complex plane, except at the negative integers. 

B. Psi digamma function 

The psi C^*) function, also known as the digamma function, is the function of a complex variable z 
obtained by differentiating the logarithm of the gamma function [20], i.e 

*(^) = ^ (3) 



It is known that ^(1) = —7, where 7 0.57 is known as the Euler-Mascheroni constant |20 



C. Pochhammer's Symbol 



The Pochhammer symbol {z)n is defined as the ratio between F(z + n) and F(z) |20|, i.e 



_ T{z + n) 
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D. Generalized hypergeometric functions 

A generalized (p, q) order hypergeometric function pFq is defined as a power series of the complex 



number z. The expression of such function can be written in the form |20| 



ai,...ap \ _ ^ (ai)n(Q2)n • • • {ap)n 



Many calculus packages such as Mathematica |21 1, Maple |22| and Matlab p3| have routines to compute 
these hypergeometric functions. Yet, theses functions have interesting properties and in some cases they 
can be converted to polylogarithms which when evaluated at the unit circle can be expressed by standard 
functions p4l. 

E. Polylogarithm functions 

For some specific values of and in ([5]), the hypergeometric functions can be converted to 



polylogarithm functions |25|. The polylogarithm, also known as the Jonquiere's function (See |26|), 



is the function defined in the complex plane over the open unit disk by 



oo ^ 
k=l 



For n = I, one simply obtain Lii{z) = — ln(l — z). The polylogarithm function is connected to the 
generalized hypergeometric functions through the relation p4| 



l,a2 ••• an+i . 

n+lFn I ;z ] Z = Lln[z) 

02 + 1,03 + 1 ... a„+i + 1 



More specifically, one has 124] 




z \ z = Li2{z), and (7) 



4F3 ^zjz = Lh{z) (8) 



2,2, 2 



F. Polylogarithm on the unit circle 

For \z\ = 1, z = e^"^^^ (0 < A < 1), the polylogarithm functions can be expressed as Li„(e''^'^^) = 
C„(A) + iSn(A), where Sn(A) = 9[Li„(eJ'^)], and C„(A) = 3?[Li„(eJ'^)]. One has the following first 
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C„ and Sn functions |[20| 



Si(A) = vr(i-A) (9) 

C2(A) = TT^d-A + A^) (10) 

D 

S3(A) = 7r3A(^-A + ^A2) (11) 



It will be shown later that (|9l)-([TT]) are the polynomial functions that appear in the expressions establishing 
the boundary condition for subharmonic oscillations occurrence in the buck converter. While the first one 
(|9]) corresponds to the peak and valley current mode control, the second and the third ones ( 10l-( 11 



appear in both voltage mode control and current mode control with dynamic compensator such as average 
current mode control 1161. 



G. The Riemann Zeta function 
It is the function defined by pO] 



oo ^ 



k=l 



For example, for z = 2, C,{2) = — . It can be proved that one has the following values of C,{2k) for even 

6 

integer arguments pO[ 



where are called Bernoulli numbers ||27|. Explicit formulas for the values C,{2k + 1) for odd integer 
arguments are not available and currently it seems that this is a very difficult problem for mathematicians. 
In fact, it is only proved that C(3) ~ 1.202 is an irrational number known as the Apery's constant. For 
engineering use, the fact that C,{k) converges quickly to 1 can be used to approximate C,{k) by 1 for 
/c > 3. 

IV. Buck converter under voltage mode controller 

The reader may now ask why all previous rather theoretical material is presented? Simply because 
these functions will appear in the mathematical expressions describing the periodic behavior and the 
subharmonic oscillations in the buck converter. The presentation of this theoretical material will allow to 
the reader a better understanding of the results presented in this paper. 
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A. Dynamic model of the power stage circuit 

Consider a buck converter under voltage mode control shown in Fig. [T]-(a) that can be equivalently 
represented by the block diagram shown in Fig. [T]-(b). For simplicity, let us suppose that = Ls + 
and Zc = rc + 1/Cs and that the load is purely resistive {R). The approach can be extended to other 
more complex models of the loads and reactive components but at the expense of more mathematical 
involvement. For the considered case, the buck regulator power stage input-to-output {5-to-v) transfer 
function G^sis) can be expressed in the following generalized form |[l| 

^ + 1 

Gvsis) = Gvso—2 — (12) 

^ + ttVs + 1 



All the parameters that appear in the expression of the transfer function ( [12] ) are given in a list of 
parameters given at the beginning of the this work. G^so is given by 

Gvso = K^eVg, where = — 

R + ri 

Theoretically, parasitic parameters can be easily ignored. However, the first order parasitic elements such 
as ESRs of the inductor and the output capacitor are included in the analysis. The effect of these parasitic 
elements may be insignificant in some cases but it may be dramatic in others. Engineering judgement 
can be done once included to examine their effects. 



B. Closed form conditions for periodic behavior and subharmonics 

Let Gc{s) be the compensator transfer function from the error Ve{t) to the control voltage Vc{t). Let D 
be the steady state duty cycle and z = e^'"^^, where A = D for trailing edge modulation and A = 1 — D 
for leading edge modulation. Let V/ be the lower value of the triangular ramp modulator signal vtri{t), 
Vm its amplitude and T its period. Let G{s) = Gys{s)Gc{s) and Gk = G{kjujs)- Let us also define 

Hk = Gk/vg. 

It can be recognized that the voltage at the input of the Z^ZcR low pass filter is essentially a square 
wave with amplitude Vg. The period of this signal depends on the dynamics of the closed loop system. 
Expanding this square wave in a Fourier series and operating on each term G^, equating the resulting 
expression to the ramp signal at the switching instants defined by the crossing between Vc{t) and vtri{t), 
conditions for different periodicity can be obtained. This approach has been applied in |28| to obtain the 
following condition for periodic behavior. 

VgiHoD + 25?[J] ^ ~.^J ' ]) - GcOVref = Vi + AVm (13) 

k=l ^ 
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where Gco = G'c(O) and stands for taking the real part. In |28|, it has been also shown that at the 
boundary of subharmonic oscillations, the following equality holds 

oo 

2vg^^{l - z'')Hk - Hj,_,] = Vm (14) 

k=l 

The system of equations ([T3|)-([T4|) can be solved for any design parameter to locate the boundary between 
the desired stable periodic behavior and subharmonic oscillations. This equation can be solved either 
graphically by plotting ( [T3| ) and (\A\ for a sufficiently high number of terms and looking at the intersection 
point, by numerical methods or analytically for some limit cases as it will be shown later. 

In 1 28 1, the series in ( 13 l-ifT?]) have been approximated by the term that involve the transfer function 
G{s) with the smallest argument or solved graphically by truncating the series. In this work, closed form 
expressions will be given for the series involved in ( [T3] ) and ([14]). To calculate the series, it is required 
to use the concept of hypergeometric series and polylogarithm functions presented in Section II. It will 
be shown later that by considering practical operating conditions, these series can be approximated by 
standard functions depending on the power stage and controller parameters and more importantly on the 
steady state value of the duty cycle D. 

1 ) Conditions for periodic behavior: From ( [T3] ), an expression for the input voltage in terms of the 
circuit parameters and the duty cycle D for the periodic regime is 

+ Vi + VmA 

^-^(^^ = H,D + e{D) ^''^ 
From (15 1 and ( [T3] ), e{D) is a small quantity depending on D given by 

siD) = (16) 

k=l 



where 9 stands for taking the imaginary part. Without £{D), ( fTS] ) represents the well-known relationship 
between the input and the output voltages of the averaged model of a buck converter. Therefore, the e{D) 
term represents a correctional factor that determines the difference between the averaged model and the 
reality of the switched system in the case of periodic regime. Note that the same approach is valid for 
determining other critical parameter values, such as the feedback gain, ramp amplitude and reference 
voltage. 

2) Conditions for subharmonic oscillations: Similarly, from the equation establishing the subharmonic 
oscillations boundary ([14]), the input voltage at this boundary is 

v,,m = ^ (17) 
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where from ([14]), S{D) is given by 

oo 

S{D) = 2^^{l-z'')Hk-Hk^,] (18) 
fc=i 

The function S{D) depends on the controller used. In the next sections, closed form expressions will be 
given for £{D) and S{D) for different controller transfer functions. These expressions will be derived 
exactly in terms of hypergeometric functions and approximately by using polylogaiithms and their special 
expressions on the unit circle using Eqs. (|9l)-([TT]). 

V. Application to a buck converter under voltage mode control 
A. Buck converter under a simple proportional controller 

It is instructive to start the application of the approach to a buck converter with a simple proportional 
controller. In this case, the controller transfer function is simply a constant for all the frequencies 
{Gc = ky). Both trailing edge modulation and leading edge modulation are studied in a unified approach. 
The results obtained for the leading edge modulation case can be applied directly to the trailing edge 
modulation case by simply changing D hy 1 — D. 

Let us consider that the output capacitor is ideal. However, the results for this simple case will be 
valid for a buck converter with a realistic output capacitor and under a single pole controller with a pole 
ojpi located at exactly the zero ojzi due to the ESR of the output capacitor. 



1) Condition for periodic oscillations: By using ^V7\ and ( [76] ), the expression of £{D) can be written 
in the following hypergeometric-function-based form 

^''''^-^[^/(l -Ed- jSr) + ^(1 +ed- jSr) 



2tt 
3 



-Ed- j£r) + ^-(1 + ed- jSr)) 



^ 2z /l,l,l-erf-je„l + erf-je, 

+ 1 -j -4^3 ; 2 I J (19) 

if £oQo ^ \ 2,2 + erf-jer,2-ed- je^ 



where Er = Eq/2Qq and Ed = a/4Qo ~ 1- ^^^^ controller, Gco = k^. Ho = k^Kg and (15l becomes 



, n^ kyVref + Vl + Va/A 

^9.1 (^) = — , n I — rrr\ — ^^^^ 
kyKiD + e{D) 

Due to the ujpi — Uzi pole-zero cancelation, this is the same expression that would be obtained for a 
single pole controlled buck converter with an ESR in the output capacitor. 
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In a practical design, ujo « ujg and ojqQq « ojg and then, eq 
becomes 



lim lim s{D) 

£o->-0 EoQo->-0 
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4F3 


^1,1,1,1^^^ 




- C(3) j 


vr 






^ 2,2,2 '^J 







and eoQo 0. Then, ([19]) 



(21) 



where C(3) is the Apery's constant ( |29|). After calculating the specific hypergeometric series, (2]_l gives 



lim lim e(-D) 



vr 



(9 [Li3(e2-^-^)]-C(3)) 



(22) 



Lis is the trilogarithm function defined in general by Q. The imaginary part of the trilogarithm evaluated 
at the unit circle is the function S3 defined in ( [TT] ) and therefore one has from (22 1 



lim lim e(-D) 



vr 



(S3(A)-C(3)) 



(23) 



For design purpose, e{D) can even be ignored because the first term in the denominator of ([20| is 
dominant. In fact it can be proved that k^K^D » e{D) for the whole operating range of the duty cycle 
D and therefore the input voltage Vg^i{D) will be related to the duty cycle D, in the case of a single pole 
controller canceling the ESR zero w^i or in the case of a simple proportional controller, by the following 
approximated expression 



^9 AD) 



kvVref + Vi + VmA 



(24) 



k^KiD 

Note that this is the same expression that one would obtain if a simple averaged model is used. For design- 
oriented analysis (|24]) is accurate enough while an approximated expression for £{D) is still available 
in terms of standard functions. Indeed, if S3 is substituted by its expression in ( [TT] ) and by using (20), 
another approximated closed form expression for Vg^i{D) is 

_ ^kyVref + Vi + Va/A) 

""^'^^ ' ^ 3kvKeD + kvKi£liT^A{l-3A + 2A^) ^ ' 

Note that both ( [24] ) and ( [25] ) are approximated expressions while the exact equation is ( [20] ). 

2) Condition for subharmonic oscillations: Following a similar procedure to that of the previous 



section, an expression for S{D) in (J^) can be obtained. For the case of single pole controller canceling 
the zero lozI due to the ESR of the output capacitor or in the case of a simple proportional controller and 
an ideal output capacitor, from (12) and ([TT), S{D) is given by the following hypergeometric-function- 
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based expression 

SAD) 



1 ^ 



1 



*(1 + ed- jsr) + ^-l- +ed- jSr)) 



-3F2 



l,l + ed-jer,l-£d-j£r 



; 2; 



(26) 



Under the practical conditions, Sq ^ and EqQo —^0, (26) becomes 



lim lim Sp{D) 



fc^eg 2C(2) + 3? 



3-^2 



1,1,1 



:z ] z 



2,2 



(27) 



By calculating the special hypergeometric series that appears in (27i, Sp{D) can be written in the 
following dilogarithm-based expression 



lim lim Sp{D) 



(28) 



Using ( fTT] ) and ( [28] ), the critical input voltage Vg^2{D) giving the boundary of subharmonic oscillation 
in terms of the duty cycle D is obtained 

Vm 



(29) 



A;,K,e2(2C(2) + K[Li2(e2-iA)]) 

Taking the real part of the dilogarithm function using ( [TO] ) and substituting A, the following expression 
is obtained for both the leading edge modulation and the trailing edge modulation modulation strategies 

Vm 1 



k^Keel TT^l - 2D{1 - D)) 



(30) 



which is the same expression obtained in 1 33 1 for the case of ideal components (re = and = 0) by 
using a very different approach. Note that ki and eo in ( |30l ) depends on the parasitic parameters. The 
critical value of the input voltage depends on both power stage and control system parameters and most 
importantly on the duty cycle D. 

From (|30|), subharmonic oscillations can be avoided in the system if the following inequality holds 



1 



D + D^ < 



Vm^s 



2vgk^K£UjQ7r'^ 



(31) 



where / is a second degree polynomial function of the duty cycle D and g is a function of the vector p of 
the system parameters. Eq. ( |3T] ) can be considered as an extended expression to ([T]) for the voltage-mode- 



controlled buck converter. It can be observed that ( 3 1 1 is invariant under the change D ^ I — D which 
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(a) Vg^i{D) and Vg^2{D) in dB (b) Exact and approximated Vg^2{D) 

Fig. 3. (a) Mesh of Vg^i{D) and Vg^2{D) in terms of the duty cycle D and the quality factor Qo where Vg^i{D) and Vg^2{D) 
are represented in dB. (b) Mesh of Vg^2{D) from l[T7} and \30\ showing that for small Qo, \30\ is not enough accurate. 



implies that the same expression is valid for both leading edge modulation and trailing edge modulation 
strategies. To validate the previous theoretical predictions, let us consider the following example. 

Example 1: Consider the well known and widely studied example of buck converter considered first 
in Q and later by other researchers in ||7j, |[8|, ||9|, p3| . The same set of parameter values will be 
considered so that the readers can make the comparison easily. Namely, L = 20 mH, C = 47 ^F, 
Vi = 3.8 V, Vm = 4.4 V, r = 400 /zs, Vref = 11-3 V and ky = 8.4. Furthermore, numerical simulations 



are not repeated here to save space. Interested readers can see 1 13 1 and 1 14 1 for both numerical simulations 
and experimental measurements. 

Figure [3]-(a) shows a mesh plot of Vg^2{D) from (17) and Vg^i{D) from (J_5l in terms of the duty cycle 
D and the quality factor Qo or equivalently the load resistance R. The intersection of the two surfaces 
is the curve of the critical input voltage at the boundary of subharmonic oscillations. In Fig. [3]-(b), the 
exact mesh plot of Vg^2{D) from (Hi is shown together with the approximated plot from ( [SO] ). From this 



figure, it can be observed that for small Qq, a discrepancy exists between the exact and the approximated 
plots. This discrepancy becomes significant for a power quality Qo approaching 1/2. For Qo < 1/2, ( |30l ) 
will give inaccurate results. In fact, the main drawback of the simplified expression based on ( |30l ) is that 
it practically does not depend on the load resistance R which is an important design parameteiQ 

'Note that the unique dependence of the simplified expression of Vg^2{D) on R is through the parameter K£ « 1. 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Duty cycle D Duty cycle D 

(c) i? = 22 n ^ Qo ~ 1 (d) i? = 50 n ^ Qo ~ 2.42 

Fig. 4. Stability boundaries of the buck converter under a simple proportional control in terms of the duty cycle for different 
values quality factor Qo or equivalently the load resistance R. Solid: exact curves. Dashed: approximated curves. 



Figure |4] shows the exact plots of Vg^2{D) and Vg^i{D) from ( [T5] ) and ( [TT] ) and the approximated plots 
from (|25ll and dSOb for different values of the load resistance R. For all the considered values of this 



parameter, the curves of Vgi{D) using ( [TS] ) and (|25]l are practically coincident. However, concerning 
Vg,2{D), it can be observed that there is a discrepancy between the exact curve and the approximated 
one for low values of the load resistance. 



The critical value of the input voltage Vg^cri is the intersection of the curves Vg^i{D) from (15 1 and 



Vg,2{D) from ( [TT] ). It can be determined by simply subtracting Vg^i{D) from Vg^2{D), equating the 
resulting expression to zero and solving for D and then substituting in Vg^i{D) or in Vg^2{D). For 
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instance, for R = 22 (Fig. |4[(c)), the critical value of the input voltage Vg^cn ~ 24.51 V (for D ^ 0.47) 
which is in perfect agreement with |[8|, p3| , |14|. For this value of the load resistance, Qo ~ 1 and the 



critical value can be also obtained accurately from the approximated expression. 



As Qq decreases, the value of intersection point Vg^cn increases as predicted in 1 13 1 by using a different 



approach based on Fillipov method and the monodromy matrix. For example, for = 5 Q., the critical 



value of the input voltage is approximately 31V which is in perfect agreement with 1 13| (See Fig. 7 in 



1 13 1). From the approximated expression of Vg^2{D) one still obtains Vg^cn ~ 24.51 V. This discrepancy 
is due to the low value of Qo- 

Note that the symmetry with respect to = 1/2 is lost for low values of the quality factor Qo which 



does not appear in (30 1. Also, the maximum value of Vg^2{D) moves to the right side of D = 1/2, when 
Qo decreases in the case of leading edge modulation while it moves to the left side in the case of trailing 
edge modulation. 

B. Buck converter under a PI controller 

A simple dynamic compensator that remove the steady state error from the system under a proportional 
control, can be obtained by adding a pole at the origin. Also a zero ujz2 is added to improve the 
phase margin. The resulting controller is a Proportional-Integral (PI) compensator and generally has 
the following transfer function 

G,{s) = k'-±^ (32) 



The added zero ujz2 is selected to be at aujQ (0.5 < a < 1.6) |31 1, |32|. Note that the total loop transfer 
function will be the same one corresponding to a Type II controller with an aditional pole ujpi canceling 
the zero ujzi due to the ESR of the output capacitor |32|. The exact expression of Vg^i{D) can also 



be expressed by ( [TS] ), where e{D) is also given by a generalized hypergeometric function and can be 
approximated by poly logarithms. However, as Gco and Hq are infinite for this type of controllers due to 
the pole at the origin, from ( [T5] ), Vg^i{D) can be approximated by Vref / D which is a very simple and 
well known expression for Vg^i{D) relating the input and the output voltages in a buck converter that 
can also be obtained from a conventional averaged model |[T|- Following a similar procedure as for the 
case of a simple proportional controller, the exact expression of Vg^2{D) in this case is also given by a 
hypergeometric function. The same expression of ( [TT] ) is valid in this case with Spi depending on a new 
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parameter = a^o and given by the following expression 

SAD) = -Ed- jer)ea - 2£oQo^(l -Ed- jSr 

2 ^AQl - 1 

+ ^'(1 + ed- jer)ea - 2Qoeo^(l + - j£r) 
+ ^(^ - - jer)ea - 2eoQo^(^ - - j£r) 
- ^(^ + £d- jer)ea + 2eoQo*(^ +ed- jSr)) 

+ + Ed - jer)£a + -Sd- j£r)£a 

+ + erf - j£r)a - - ffd - jer)a 

2jz{j + ea) (l,l,'2-jea,l-ed-j£r,l + ed- ^ . 
+ -I 1 , . — 5-f^4 I J (33) 

^ iJ+-?eo(9o \ 2,l-ie„,2-erf-jer,2 + erf-je^ ' 




Note that when a — 0, — 0, and Spi{D) — 5p(D). To obtain a more simple expression for Vg^2{D) 
in the case of a type II and PI controllers, let us consider eq — ^ and EqQq — )• 0, which is the case in 



almost all practical designs. Therefore, in terms of poly logarithms, ( |33l ) becomes as follows 

lim lim SAD) = KKiel{^[U2{e^^'^)] + a'i))-ael%U^{e'^^^)] (34) 

It can be noted that, as >> (Jq (eo — ^ 0), the contribution due to the trilogarithm term is actually 
negligible. The real part of the dilogarithm and the imaginary part of the trilogarithm have been already 
obtained before. Therefore, for the case of Type 11 and the PI controllers with a trailing edge modulation 
(A = D), the approximated critical value of the feedback gain has the following expression 

' ^ KtVgel^^^{l-2D{l-D))+aea^^D{l-?,D + 2D^) ^ ' 

Similar analysis can be applied to obtain the expression corresponding to a leading edge modulation 
strategy. Note that the symmetry with respect to D = 1/2 is lost although this cannot be appreciated 



since Eq is very small in practice. Equation ([35]) can be considered as an extension of ( |30| ) for the case of 
the Type II and the PI controllers. The next example validates the results for a buck converter for these 
two different controllers. 

Example 2: Consider a buck converter with voltage mode control and trailing edge modulation studied 
in p3) whose value of parameters are oriented to miniaturization. The considered fixed parameters in 
this example are: = 1 17, L = 30 nH, C = 50 nF, fs = 50 MHz, Vi = Q, Vm = I ^, Vg = 2, N and 
a = 1/2. Two cases will be studied. In the first one, ideal output capacitor and inductor are considered 
and the system will be under a PI control. In the second case, parasitic elements will be included. While 
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for the ESR of the inductor, only static performances are changed, in the case of an ESR Tc of the 
output capacitor, a zero ujzi is introduced in the dynamic model of the system. This zero is canceled 
with an additional ujpi of the type II controller. Note that although the dynamic effect of the ESR output 
capacitor zero is canceled, its static effect still appears in the expression of Vg^2{D) since eq depends on 
this parameter. 




Duty cycle D Duty cycle D 



(a) Tc = 0, r£ = 0, a = 1/2 (b) = 0.05 Q,, n = 0.1 a = 1/2 

Fig. 5. Stability curves of the buck converter of Example 2. (a) fc.u in terms of the duty cycle D for an ideal converter 
fci>,cri ~ 9.7. (b) fc„ in terms of the duty cycle D taking into account parasitic elements. fci,.cri ~ 10.7. 



1} Ideal components: The system in this case is under a PI controller. Furthermore, in order to do 
not induce slow scale instability, let us select a such that the zero a;^2 = ^^0/2 (a = 1/2) of the PI 
controller be below ujqQq |33J. Therefore, it will be below ojg because ujg > (^oQo in a practical design. 
Figure [5]-a shows the stability curve ([35]l in terms of the proportional gain and the duty cycle D. From 



p5} , the critical value of the feedback gain for the previous set of parameter values with Vref = 1.5 V 
(D = 1/2), is ky^cvi ~ 9.7. 

2) Taking into account parasitic elements: Let us consider that rc = 0.05il and = 1 mD.. The system 
in this case is under a Type II controller. Its additional pole Wpi is selected at exactly ujzi = l/vcC. As 
before the zero W22 = a;o/2 of the PI controller is below ujqQq. Figure [5]-(b) shows the stability curve 
in terms of the proportional gain ky and the duty cycle D from ( [35] ). The critical value of the feedback 
gain for the previous set of parameter values with Vref = 1-5 V (L> = 1/2) is k^^cri ~ 10.7. 

Figure [6]-(a)-(d) shows the time domain numerical simulation for the two previous cases for value 
of the feedback gain just before and just after subharmonic oscillations occurrence. These numerical 
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(a) rc = Q,ri = Q,a = 1/2, fc„ = 9.5 



(b) r-c = 0.05 VL,ri = l mO., a = 1/2, fc„ = 10 
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(c) rc = 0, rf = 0, a = 1/2, fc„ = 10, 



(d) rc = 0.05 n,ri = l mO,, a = 1/2, fc„ = 10.8 



Fig. 6. Dynamic behavior of the buck converter of Example 2 just before and just after subharmonic oscillations occurrence, (a) 
periodic behavior of the system with ideal components under a PI controller, (b) Subharmonic oscillation of the system under a 
PI controller with ideal components, (c) Periodic behavior of the system under a Type II controller taking into account parasitic 
parameters, (d) Subharmonic oscillation of the system taking into account parasitic parameters under a Type II controller. 



simulations are in perfect agreement with the theoretical derivations. 

VI. Design-oriented prediction of subharmonic oscillations 

A. Crossover frequency as a Figure of Merit 

The closed form results presented previously can be interpreted in terms of several Figures of Merit like 
for example the ripple amplitude in |33|. The results can be also interpreted by the crossover frequency 
and the associated phase margin which depends on the controller used. The crossover frequency ujc is 
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defined as the frequency where the modulus of the loop gain crosses dB {\T^{juJc)\ = 1)- For instance, 
for the buck converter under the proportional, PI and Type II controllers, the crossover frequency can be 
approximated by 

u.^uo^— (36) 

Note that increasing the proportional gain ky or decreasing the modulator amplitude Vm will raise the 
crossover frequency but this can also cause subharmonic oscillations. From ([30]) and ( |36l ), the boundary 
conditions for this phenomenon, in terms of the crossover frequency can be written in the following 
general form 

« u^,p{D) (37) 

where p{D) depends on the controller used. For the case of proportional, PI and Type II controllers and 
if a;o << ujs and ooqQo << u^s, from ([30]), p{D) can be approximated by 



From ( [36] ), the condition for the system to be free form subharmonic oscillations will be 

Uc < 0Js\^p(D) ■■= oJcMD) (39) 



In 1 30 1, the upper limit of the crossover frequency is discussed. For an ideal converter it was recommended, 



without demonstration, that an upper limit is w^/S. Eq. (39 1 establish a theoretical upper limit for avoiding 
subharmonic oscillations in the buck converter. 

For instance, in Example 1, taking Vg = 23.5 V and R = 22 Q,, the value of the duty cycle is 
D » 0.5 and from ([39]), ujc = 6.95 krad/s < a;c,cri(^) = 7.069 krad/s and the system is 
stable. For Vg = 25.5 V and R = 22 Q, the value of the duty cycle is D ^ 0.48 and from ( |36j ), 
cjc = 7.24 krad/s > LOc,cvi{D) = 7.059 krad/s and therefore the system exhibits subharmonic 
oscillation for this set of parameter values. 



B. Phase margin as a Figure of Merit 

Traditionally, the phase margin ipm is used, as a test on the loop gain T-u{s), to determine the stability 
of the closed loop system Tci(s) = T„(s)/(1 + Ty{s)) in the averaging context. The phase margin is 
determined from the phase of Ty{s) at the crossover frequency lOc , by 

(Pm = 'n- + Z.Ty{jujc) (40) 
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If |Tl,(s)| crosses dB only once, then the closed loop system T^(s)/(1 + Ty{s)), the system Tci{s) 
will not present slow scale instability that can be predicted by a simple averaged model. Traditionally 
a phase margin of at least 45° is required to guarantee an acceptable system response |[T|, | [3T| , | [32} . 
Furthermore, to avoid subharmonic oscillations, the phase margin must be greater than Wc,cri(-D) in 



(39l 



Note that the phase margin will depend on the type of the controller used and the system parameters. 



From p9] ) and ( |40| ), the general expression of the critical phase margin can be written in the following 
form 

(Z?) «7r + Zr,(i^,,criP)) (41) 

Subharmonic oscillations can be avoided if tpm < ^m,cn{D). For (pm.cn{D) < 45°, a well designed 
converter based on the loop gain analysis of the averaged model will not exhibit subharmonic oscillations. 
For ipm,cTi{D) > 45°, subharmonic oscillation may still occur even for a well designed converter. 

Conclusions 

Subharmonic oscillations in peak-current-mode-controlled converters are well documented by different 
analytical methods. However, in voltage mode control, this phenomenon has been only characterized by 
using numerical simulations or mainly based on abstract mathematical analysis using a discrete time 
model or the Fillipov method. Both approach do not allow to relate the results to concepts that are 
familiar to power electronics engineers. Explicit expressions for conditions of subharmonic oscillations 
occurrence have been unavailable for many years. 

It is widely known that for peak current mode control without ramp compensation, for example, the 
stability condition is D < 1/2. This expression is shown to be a special case of the general results 
presented in this paper. In general, it can be conjectured that the stability condition is in the form 
f{D) < g{p), where / is a function of the duty cycle D that can be approximated by a polynomial 
function and g( is a function of the vector of system parameters p. It was also shown that the structure of 
the loop gain offers some additional insights on issues related to the stability of the system. Interestingly, 
the degree of the polynomial function / is extremely related to the relative degree of the total loop 
transfer function. For instance, for peak current mode control, the relative degree of the total loop gain is 
1 and the degree of / is also 1. For voltage mode control, the relative degree of the loop gain depends on 
the compensator. For proportional controller, single pole controller canceling the ESR zero of the power 
stage and for PI controller, the relative degree of the total loop is 2 and / can also be approximated 
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by second degree polynomial function. The same can be conjectured for total loop transfer functions 
corresponding to other control schemes such as average current mode. 

Based on the approach presented in this study, critical values of the system parameters can be located 
accurately. Only in special limiting cases, an analytical solution can be obtained. Analytical methods end 
at a certain point and have to be succeeded by numerics in the general case. However, hypergeometric- 
based functions have some physical insights on how the system parameters affect the dynamics and how 
approximations should be done to predict subharmonic oscillations in an analytical way in the case of a 
quality factor Qo > 1/2. Furthermore, by guaranteeing the stability for a high Qq case, the stability is 
also guaranteed for the low Qo case. 

After explicitly obtaining the conditions for subharmonic oscillation, the results have been reformulated 
in terms of Figures of Merit widely used by the power electronics community, such as the crossover 
frequency, and phase margin. The critical expressions for these Figures of Merit can be obtained in 
terms of converter parameters and its controller. The ability to predict subharmonic instability using pure 
algebraic equations directly expressible in terms of such Figures of Merit is a major advantage of the 
proposed approach. 

Works are in progress to extend the approach to other switching converter topologies and other control 
schemes. The results will be reported in a further study. 
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